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ABSTRACT 



Anomalous flux ratios between quasar images are suspected to be caused by substructures in lens galaxies. We present new deep and 
high resolution H and Ks imaging of the strongly lensed quasar SDSS J0924+0219 obtained using the ESO VLT with adaptive optics 
and the Laser Guide Star system. SDSS J0924+0219 is particularly interesting as the observed flux ratio between the quasar images 
vastly disagree with the predictions from smooth mass models. With our adaptive optics observations we find a luminous object. 
Object L, located ~ 0.3" to the North of the lens galaxy, but we show that it can not be responsible for the anomalous flux ratios. 
Object L as well as a luminous extension of the lens galaxy to the South are seen in the archival HST/ACS image in the F814W filter. 
This suggests that Object L is part of a bar in the lens galaxy, as also supported by the presence of a significant disk component in 
the light profile of the lens galaxy. Finally, we do not find evidence for any other luminous substructure that may explain the quasar 
images flux ratios. However, owe to the persistence of the flux ratio anomaly over time (~ 7 years) a combination of microlensing and 
milli-lensing is the favorite explanation for the observations. 

Key words. Gravitational lensing: strong. Galaxies: quasars:individual:SDSS J092455. 87+021924.9 



1. Introduction 

Strong gravitational lensing is now part of the standard 
toolbox of the astrophysicist. It co nsists in a p ower- 
ful test for cosmological models (e.g. Barnabe et al.| [ 
Suvuetal. 2010; Coe & Moustakas 2009; FedeH et alJ 



Metcalf & Zh"aoll2002l;lKoopmans et al .1120021; iMao & Schneideii 



1998). 



2011 



2008 



Ipedeli & BartelmannI 120071) and allows us to study the distri- 
bution of luminous and dark matter in galaxies as well as 
the evolution of their properties with redshift (e.g. Ruff et all 
12011'; 'Faureet al."2011t iTortora et al.ll2010t lAuger et al.ll2'oT(if 
iKoo pinans e t al. 2006). 

Such studies require a model for the total projected potential 
well of lensing galaxies, which is often described as a smooth 
2D elliptical profile, with or without contribution of intervening 
objects along the line of sight. However, the increasing accu- 
racy of the observational constraints provided by high resolu- 
tion imaging, has quickly shown the limitations of such simple 
models: in quadruply imaged quasars, smooth models are often 
unable to account simultaneously for the milli-arcsec astrome- 
try a nd for the near- and mid-IR flux ratios of the quasar im- 



ages ('X^iettietaD "20103; 'Chantry et al.'l 20 Ot lYoo et al 
[Biggs et al. 2004; Kochanek & Dalai 2004t ' 



Keeton et al. 



2005 



2003 



The reason for the observed discrepancies between the 
model predictions and the measurements may be the pres- 
ence of substructures in the halo of lensing galaxies. Such 
small deviations to the smooth potential might affect signifi- 
cantly the predicted flux ratios while barely changing the as- 
trometry of the quasar images. This is especially true when 
a quasar image is a saddle poi nt of the arrival time surface 
jSchechter & Wambsganssll2002h . Strong lensing therefore of- 
fers a sensitive way to indirectly detect and weight substruc- 
tures in galaxy halos, whether luminous or not (e.g. Vegetti et aP 
l20T0bi: ISuvu & Halkolall2010t iMcKean et al]l2007h . 

In the present work, we have searched for substructures 
in the halo of the lens galaxy towards the quadruply imaged 
quasar with t he mo st anomalous flux ratios, SDSS J0924+0219 
dlnada et aTl l2003h . at Zquasar — 1-524, discovered in the 
Sloan Digital Sky Survey (SDSS, lYork et al. 2000). Hubble 
Space Telescope (HST) spectra of the system rev eal tha t one 
of the quasar images, labelled D in ilnada et al.l (l2003h (see 

al.1 120061) ^1 



* Based on observations obtained with the ESO VLT at Paranal ob- 
servatory (Prog ID 084.A-0762(A); PI: Meylan). Also based in part on 
observations made with the NASA/ESA Hubble Space Telescope, ob- 
tained at the Space Telescope Science Institute, which is operated by 
the Association of Universities for Research in Astronomy, Inc., under 
NASA contract NAS 5-26555. These observations are associated with 
the CASTLES (Cfa-Arizona Space Telescope LEns Survey) survey (ID: 
9744, PI: C. S. Kochanek) 



also Fig. |2]l, is extremely faint (iKeeton et al.l 120061) . It is also 
a saddle point in the arrival time surface, hence making it 
mo re lik ely to be demagnified owe to micro- and milli-lensing 
(ISchechter & Wambsganssl l2002). The microlensing hypothesis 
provides a satisfying explanatio n for the observed flux ratio 
anonialy ( Bate et al. 201 1; Medi avilla et al.ll2009l; iMorgan et"al] 
120061 : iKeeton et al.ll2006) . However, this flux anomaly lasts for 
at least 7 years, which is significantly more than the expected du- 
ration of microlensing fluctuations in this lens (0.39 years), and 
larger than the duration averaged over a sample of 87 lenses, i.e. 
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7.3 month, (iMosquera & Kochanekll20TTI) . It is therefore con- 
ceivable that not only stellar-microlensing takes place in this 
system but also that a more extended/massive substructure also 
contributes to the observed anomaly. A deep spectrum of the 
lens galaxy was obtained with th e ESQ Very Large Telescope 
(VLT) and the FORSl instrument jEigenbrod et alj |2006). lead- 
ing to ziens = 0.394 + 0.001 from stellar absorption lines. In addi- 
tion Eigenbrod et al. (2006) notice numerous elongated features 
within and around the Einstein ring in the HST images. This ob- 
servation leads them to conclude that the lensed source may be 
double. 

In order to investigate further the nature of these features 
and to look for possible faint and small satellites to the lensing 
galaxy, we conduct deep Adaptive Optics (AO) imaging observa- 
tions of SDSS J0924+0219 with the VLT. The paper is organized 
as follow: in Section|2l we present the VLT AO observations and 
the HST dataset. In Section [3] we detail our findings of a lumi- 
nous substructure in the AO images. In Section]!] we show that 
the substructure is also present in the optical, from HST/ ACS 
observations. In Section]5] we discuss the possible nature of the 
substructure. Finally, we present our conclusions in Section]6] 

Throughout this paper, the WMAP5 ACDM cosmology is 
assumed (Q.^ = 0.258, Qa = 0.742, Hq = 72 km s"' Mpc"'). 
All magnitudes are in the AB system. 

2. The datasets 

2. 1 . Adaptive optics imaging witli tiie VLT and NACO 

We observed SDSS J0924-I-0219 with the near-infrared camera 
CONICA, which is mounted on the Adaptive Optics system 
NAOS, NACO for short, installed at the Nasmyth B focus of 
the VLT-UT4 at the ESO Paranal Observatory, Chile. Our ob- 
servations were obtained on 2010 March 12, and 14 (H- and 
Ks-bands) and on 2009 December 14 (Ks-band), using the laser 
guide star facility. The star U0900-06393922 (V = 16.5mag), 
located 47" away from our target, was used to correct for "tip- 
tilt". The S27 CONICA camera has a pixel size of x = 27.053 + 
0.019 mas and a field-of-view of 28" on a side. A total of 30 
(41) exposures of 3 x 60s (NDITxDIT) was obtained in the Ks- 
(H-) band, with airmasses in the range sec(z) = 1.1 - 1.3. The 
weather conditions were photometric. 

The reduction procedure follows the steps exposed in 
ISluse et al. (2008), i.e., we first subtract dark frames to all the 
near-IR science images and we apply standard flat-fielding. We 
use twilight flat-fields obtained at most 4 days before the obser- 
vations. We check for the temporal stability of these calibrations: 
flat-fields obtained 4 days apart are undistinguishable while flat- 
fields obtained 4 months apart show differences < 2%. 

The frames are divided in subsets of continuous obser- 
vations. Each frame is sky-subtracted using the the xdimsura 
IRAlQpack age. Each image is visually inspected to exclude data 
with unreliable sky subtraction, with residual cosmic -rays, or 
with too heavily distorted Point Spread Function (PSF). Finally, 
we construct two deep images by combining 27 (28) frames in 
the Ks-band (H-band) using the xnregistar task. The total 
equivalent exposure time is 4860 s in the Ks filter and 5040 s 
in the H filter. 

The PSFs of the combined images has a Full Width at Half 
Maximum of FWHM ~ 0.25" in the H-band and FWHM ~ 

' IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 



Fig. 1. PSFs created for the NACO Ks-band image. Left: zero 
order estimate of the PSF, i.e., PSFq (see text). Right: PSF5, dis- 
played here with the same grey scale as PSFq (logarithmic scale). 
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Fig. 2. NACO Ks-band image of SDSS J0924-H0219. Left: com- 
bined frame, whe re the quasar images are labelled following 
llnada et al same image where the quasar images 

have been fitted using PSF5 and subtracted using Galfit. The 
residuals at the quasar positions are negligible. 



0.20" in the Ks-band, with a moderate ellipticity, e ~ 0.1. 
Based on the observations of standard stars on 2009 December 
13 and on 2010 March 13, we derive photometric zeropoints: 
ZP(Ks) = 23.05 + 0.15 mag s"' and ZP(H) = 24.06 + 0.10 mag 
s '. The large uncertainties on the zeropoints reflect mainly 
the temporal variations of the low fre quency wings of the 
PSF and angular anisoplanatism (see e.g. Esslinger & Edmunds! 
1998). Extinction corrections of 0.034 mag in the H-band and 
0.043 mag in the Ks-band are used to flux-calibrate the science 
frames (iLombardi et al.l]20Tlh . 

2.2. HST optical imaging witli fhe ACS 

Optical HST images of SDSS J0924-H0219 are available from 
the archives (ID: 9744, PI: C. S. Kochanek). These data con- 
sist of F814W and F555W drizzled images obtained with the 
Wide Field aperture of the Advanced Camera for Surveys (ACS). 
The four F814W drizzled images were initially taken on 2003 
November 18 and 19 with an exposure time of 574 s each. The 
four F555W frames were observed in 2003 November 18 with 
a total exposure time of 547 s each. The images have been only 
recently drizzled (April 2010, iKoekemoer et al.l20p2|). The non- 
drizzled dataset was presented for the first time bv^Keeto n et al.l 
(2006) and used in Eigenbrod et al. (2006). In the following, we 
use the latest drizzled images directly from the archives. The 
pixel size in these images is 0.05". 
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3. Analysis of the NACO images 

Our goal with the AO images is to detect any possible faint 
and small satellite to the lens galaxy. Given the spatial scales 
involved, such satellites are hidden in the glare of the quasar im- 
ages or of the lens galaxy itself. It is therefore mandatory to esti- 
mate properly the AO PSF and to subtract the light contribution 
of the quasar images and of the lens galaxy from the data. 

3.1. PSF construction 

With the small field of view of the NACO instrument, no bright 
star is available to measure the PSF. We therefore create a PSF 
using the quasar images themselves. For that purpose, we use the 
following iterative method. 

First, we estimate a guess PSF by combining the images of 
the 3 brightest quasar images. Lets call it PSFq. The 3 quasar 
images are aligned and normalized to their peak intensity. PSFo 
is built by assigning, at every pixel location, the lowest pixel 
value amorig the 3 quasa r images. PSFq is then used in Galfit 
dPeng et all l2010i l2002l) to fit and subtract the 4 quasar im- 
ages. The resulting image obviously shows strong residuals at 
each quasar position. We median-combine these residuals, pre- 
viously normalized and aligned. The result of this operation is 
then weighted by the square root of PSFq, and added to the latter 
to build PSFi, our next estimate of the PSF. This relatively arbi- 
trary weighting allows us to minimize the noise contamination 
of the pixels located far away from the PSF center, i.e., where 
PSFo is already a good approximation of the true PSF. 

Second, we fit the galaxy light profile using a Sersic plus 
an exponential disk, in addition to the 4 quasar images. This 
light model is convolved by PSFi during the fit. The fit of the 
lens galaxy is not possible in the first step of the PSF creation 
because the poor quality of PSFo prevents Galfit to converge 
properly. We note that the use of an exponential disk in addition 
to the Sersic profile to model the lens galaxy is crucial to obtain 
acceptable fits. The residual image obtained with the improved 
PSF and lens model is used to compute a new estimate of the 
PSF, that we call PSF2. 

We then repeat this second step until the residuals to the fit 
at the quasar positions can not be improved anymore. A typical 
number of 5 cycles is necessary to achieve this goal. We show 
in Fig. [U the PSF obtained after 5 iterations, for the Ks-band 
image. As mentioned in Section 12.11 the PSF ellipticity is e ~ 
0.1. The apparent elongation of PSF5 in Fig.[T]is noticeable only 
at large radius and is visible because of the logarithmic scale. 
It is also typical of the NACO PSF shape, as observed in other 
targets of our observational program. Fig. |2] shows the Ks-band 
image after subtraction of the quasar images using PSF5. We do 
not observe systematic residuals under the quasar images, hence 
suggesting that PSF5 is reliable. We do not detect any luminous 
object under the quasar wings. 

Note that the 4 quasar images are located at very different po- 
sitions on the Einstein ring. The quasar host galaxy has therefore 
very different Position Angles (PA) at the four quasar positions. 
This fortunate configuration has the consequence to minimize 
the quasar host galaxy contribution when median-averaging the 
residuals at each step of the PSF construction. 

3.2. Quasar image fits 

The quasar image fits include light from the ring (Fig. |2]l. This 
prevents us from measuring accurate photometry for the quasar 
images. Rough estimates of the magnitude differences between 
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Table 1. Relative position of the quasar images, of the lens 
galaxy and of Object L to image A, as measured in the 
ACS/F814W deconvolved images. The quasar images are la- 
belled as in Fig. |2l The astrometric uncertainties for the quasar 
images reflect the internal error of the deconvolution method. 
For the lens galaxy they represent the dispersion in the Sersic 
centroids for different acceptable fits. 





ARA 


ADEC 


A 


0.0 ± 0.001 


0.0 + 0.001 


B 


+0.068 + 0.001 


-1.804 + 0.001 


C 


-0.957 ± 0.001 


-0.701 + 0.001 


D 


+0.555 ± 0.006 


-0.422 + 0.006 


Lens 


-0.176 + 0.020 


-0.860 + 0.020 


Obj.L 


-0.15 + 0.02 


-0.52 + 0.01 



images A and D are A(Md - Ma)h ~ 3.3 mag and A(Md - 
Ma)k, ~ 2 -2 mag. They confirm the flux anomaly observed by 
llnada et aP (l2003i) . even at longer wavelength. 

3.3. Gaiaxy liglit prof He fit and subtraction 

We use Galfit to fit and subtract the lens galaxy from the data. 
We are primarily interested in a clean subtraction of the lens 
galaxy light within the Einstein radius. To do so, we mask the 
area of the Einstein ring (quasars and quasar host galaxy). The 
use of masks ensure that Galfit does not try to fit the ring as 
a component of the lens galaxy. However, it can be a source of 
inaccuracies, for example if the zone masked encompasses the 
galaxy eff'ective radius. 

The PSFs obtained in Section 13.11 are used to convolve the 
analytic profiles in Galfit in the H and Ks bands. The best fits 
are obtained when the lens galaxy is described as the sum of 
two Sersic profiles (see Fig.O. The first one has a very shallow 
index, ni - 0.13 + 0.04 and an effective radius Retf.i - 0.90 + 
0.01". The second profile is much steeper, with n2 - 7.17 + 
0.65 and has a smaller effective radius, Retf,2 = 0.35 + 0.04". 
In the fits, the two profiles are separated by 61^2 = 30 + 8 mas 
and have PAi_2 ~ -80° (positive North to East). Their ellipticity 
is very small, ei 2 ^ 0.1. The values given above are averaged 
over the two near-IR bands and the error bars correspond to the 
dispersion between those values. The agreement between the fits 
in the two band is remarkable. The total magnitude of the lens 
galaxy is Ks = 15.75 + 0.15 mag and H = 16.24 + 0.10 mag. 

According to our image decomposition with Galfit the 
lens galaxy contains a significant disk-like structure, as sug- 
gested by the necessity to include a shallow Sersic profile in 
the fit. This is also sug gested by earlier results in 3 HST bands 
(Eigenbrod e t al.ll2006h . We will discuss further this finding in 
Section|5] After subtracting our model of the lens galaxy (Fig.O, 
we detect a faint but significant object North to the galaxy's cen- 
troid. The detection is significant both in the H and Ks filters : 
Ks = 22.20 + 0.15 mag, H = 23.23 + 0.10 mag, which translates 
in Lv ~ 4 ± 1 X 10'' Lq,v. This luminosity is K-corrected assum- 
ing an elliptical spectral energy distribution. We also corrected 
for extinction and evolution. In the following, we refer to this 
possible substructure of the lens as Object L. 
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Fig. 3. VLT/NACO H- and Ks-band images of SDSS J0924+0219 after subtraction of the lens galaxy with Galfit. North is to the 
top. East is to the left. The crosses indicate the fitted position of the lens. In both bands we identify a luminous residual not fitted by 
the Sersic profiles (Object L) located about 0.27" to the North of the lens center. 



ACS/F814W (O.OS'Vpix) 



Deconvolved image (0.0125"/pix) 




#1 



#2 



#3 



" Object L 



#4 



Fig. 4. HST/ACS image of SDSS J0924+0219 in the F814W filter Left: original drizzled image obtained from the archive. Middle: 
deconvolved image. The pixel size is 4 times smaller than in the left panel. Right: deconvolved image, where the lensing galaxy 
profile has been fitted and removed using Galfit. The arc-like features #1-4 are discussed in the text. North is to the top. East is to 
the left. 



4. Analysis of the HST images 

The luminous Object L discovered in the NACO images is also 
visible in the ACS/F814W drizzled images. As the field of view 
of the ACS is much wider than that of NACO several stars are 
available to build a reliable PSF and to spatially deconvolve 
the i mages. To do so, we use the MCS algorithm ( Maeain et alJ 
11998 ). The required PSF is built using 3 stars located within a 
2' radius from the lens. The spatial resolution after deconvolu- 
tion is FWHM=0.025". The adopted pixel size in the decon- 
volved data is a fourth of the original data, i.e, 0.0125". The 
original F814W drizzled image and its deconvolved version are 



displayed in Fig. |4] We do not consider here the much shallower 
F555W images. 

We then fit the deconvolved light profile of the lens galaxy, 
using a mask at the position of the ring and of the quasar images. 
As was the case with the VLT/NACO data, two Sersic profiles 
are necessary to obtain good results. The parameters of the best 
fit Sersics are : nj = 0.06, Rett i = 0.43", PAj = 9° for the 
shallow profile, and nj = 0.58,' Reff,2 = 0.14", PA2 = -20° 
for the steeper profile. Their ellipticity is low (ei 2 S 0.1) and 
their centroid are separated by <5i_2 - 4 mas. The total lumi- 
nosity of the lens galaxy is F814W - 19.36 + 0.01 mag. The 
profiles are very different from those fitted in the NACO images. 
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This might owe to the different depths, wavelength ranges and 
stellar populations probed by the data, but also to the different 
masks used during the fits and to the fact that the HST image 
(contrary to the NACO image), was preliminarily deconvolved 
before the use of Gal fit. As a consequence the photometry 
of the lens galaxy is possibly affected by systematics errors of 
several tenth of a magnitude. This result illustrates well the dif- 
ficulty to measure the light profile of a lens galaxy. However, 
we notice that both the HST and NACO data require a disk- 
like component in the fit.In addition to Object L, clearly visi- 
ble in the HST optical image, we show four arc-like features in 
Fig- El Features #1, #2 and #3 have been discovered earlier in the 
non-drizzled image (Eigenbrod et al. 2006). They were success- 
fully modeled as a secondary set of images at the quasar redshift. 
Feature #4 is identified in the ACS/F8 14W drizzled image for the 
first time. The astrometric measurements for the quasar images, 
the lens galaxy and Object L are given in Table [T] We measure 
F8 14W = 25.50 + 0.02 mag for Object L, which translates into a 
rest frame luminosity Ly ~ 9.4 x 10^ Lq v (see also Section l33] l. 

Finally, we note that the orientation of the HST image is 
about -70°. Therefore, Object L can not be a spike produced, 
e.g., by quasar image A or by the sharp nucleus of the lens 
galaxy. Removal of these spikes is in fact the main motivation 
for carrying out image deconvolution. In addition. Object L is 
present both in the space image and in the NACO images, i.e., in 
3 independent bands with two very different instruments, mak- 
ing it unUkely to be an artifact in all three datasets. 

5. The nature of the luminous substructure 

The faint Object L seen in the VLT and HST images may be 
interpreted in several ways. It can be (i) a satellite of the lens 
galaxy, (ii) the fifth central image of the lensed quasar, (iii) the 
lensed image of a source unrelated to the quasar, or (iv) the 
brightest part of a bar in the lens galaxy. In the following, we 
test the different hypothesis using simple lens modeling. 

(i) Our model consists in a smooth analytical mass profile to 
which we add a small perturber using the Lenstool software 
(Kneibetal. 1993; Julio et al. 2007). The main component of 
the model is a Singular Isothermal Ellipsoid (SIE) with an ad- 
ditional contribution of an external shear. The parameters to be 
fitted are the SIE position, its velocity dispersion, and the exter- 
nal shear strength and direction. The observational constraints 
are the positions of the quasar images available from the HST 
images (Table [1]). In addition, we use the ellipticity measured 
as a prior for the ellipticity of our mass model, i.e., we choose 
the SIE ellipticity to be lower than 0.1. This also implies that 
the model depends only weakly on the exact value of the el- 
lipticity and PA. A flat prior is also given on the SIE position, 
using the astrometric error bars of Table [T] No priors are used 
for the other parameters. The best fit leads to - 02 and 
is obtained for a SIE with crie„s = 215 km s~', which converts 
into a total mass Miens = 1.55 x 10" Mq. The external shear is 
moderate, y = 0.02, with an position angle 6y = 78°. A simple 
smooth model therefore accounts well for the astrometry of the 
quasar images, as also found by previous studies. However, the 
predicted flux ratio between the quasar images (not included in 
the;^^^ calculation) are not well reproduced at all. 

We then add a perturber to the overall mass potential. We 
describe it as a Singular Isothermal Sphere (SIS) at the position 
of Object L. This adds extra free parameters to the models. In 
order to keep at least one degree of freedom in the fit, we suc- 
cessively fix the position of the lens galaxy to 25 different values 
within its astrometric uncertainty, and repeat the modelling. The 



remaining free parameters are the SIE and SIS velocity disper- 
sions, the SIS position, and the strength of the external shear, 
i.e., we have 5 free parameters vs. 6 observational constraints 
provided by the astrometry of the quasar images. 

We then adopt flat priors on the SIS position using the as- 
trometric error bars of Object L given in Table [1] The best fit 
has;\f^ = 5.4. The corresponding velocity dispersion for Object 
L is cr£ - 0.3;^Q3 km s"'(l-cr error bar), which translates into 
a mass of Ml < 2 x 10^ Mq (1-cr). The velocity dispersion of 
the main lens is unchanged and the external shear is very little 
affected by the perturber, y = 0.018. The mass we find for the 
perturber translates in a total mass-to-light ratio lower than unity 
in all bands. If we take the 3-cr error bar on the velocity disper- 
sion of the SIS as an upper limit the mass-to-light ratio of Object 
L is barely equal to unity. It is therefore unlikely that Object L 
is a satellite galaxy or even a globular cluster within the halo of 
the lens galaxy. 

In addition, we note that introducing a perturber in the 
model at the position of Object L does not allow to reconcile 
the astrometry of the quasar images and the observed flux ra- 
tio b etween quasar im ages A, D and C as given in Table 2 in 
Eig enbrod et al.l (l2006l) . Adding more mass in Object L degrades 
the;i'^ of the fit, i.e., it affects the astrometry, and still does not 
allow to reproduced A/D flux ratio. 

(ii) An alternative interpretation is that Object L is the 5* 
central image of the quasar, like in IWinn et al. (2004). Such a 
central image is produced when the mass profile of the lens 
galaxy is very shallow or truncated in its center We describe 
this feature by introducing a core radius in our SIE model but 
we fail in predicting a central image at the position of Object 
L whatever is the size of the core radius. The hypothesis of a 
central image can therefore safely be ruled out. 

(iii) If Object L is the image of a source at the redshift of 
the quasar, we should see its counter image somewhere south to 
quasar image B and this counter image should in principle be 
brighter than Object L. We do not see any evidence for a counter 
image in the PSF-subtracted AO data nor in the deconvolved 
HST data. Object L is most probably at a redshift different from 
that of the quasar If it is at a lower redshift than the quasar, we 
might see it in absorption in the quasar spectrum. None of the 
spectra available to us (i.e., SDSSQ and VLT from Eigenbrod et 
al. 2006) show any trace of absorption lines at a different red- 
shift than that of the lens galaxy. In addition, the very deep VLT 
spectrum of the lens galaxy includes Object L in the slit and does 
not show any sign of objects at multiple redshifts. We therefore 
conclude that there is no evidence for Object L to be at a redshift 
different than that of the lens galaxy and that there is no evidence 
for Object L to be lensed. 

(iv) The last simple explanation is that Object L is part of the 
lens galaxy but that it does not affect much the lensing models. 
In the deconvolved ACS image, we detect a faint feature almost 
aligned with Object L and the core of the lens galaxy. This fea- 
ture, labelled "#4" in Fig.H] and Object L may reflect the pres- 
ence of a bar in the lens galaxy. This is also supported by the fact 
that we need a disk component to model the light profile of the 
lens galaxy, both in the optical and in the near-IR. 



6. Conclusions 

Motivated by the fact that lensed quasars with anomalous flux 
ratios can be the signature of massive substructures in the 
lens galaxy, we have carried out deep VLT/AO imaging of 
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SDSS J0924+0219, the quadruply imaged quasar with the most 
anomalous flux ratios known to date. 

With a total of 1.35 hours of exposure in the H-band and 
1.4 hours in the Ks-band using the VLT NACO system and the 
Laser Guide Star facility, we discover a luminous object ~ 0.3" 
to the North of the lens galaxy. This object, that we call Object 
L, is also seen on archival HST/ACS images in the F814W filter. 

We investigate the nature Object L using a smooth SIE 
model for the lens galaxy with an additional SIS perturber at 
the same redshift and located the position of Object L. We show 
that such a model can not explain the anomalous flux ratios in 
SDSS J0924+0219 without also affecting the astrometi-y of the 
quasar images. The parameters found for Object L that allow 
a good fit of the astrometry all yield total mass-to-light ratios 
smaller or barely equal to unity. 

We find evidence in the HST ACS/F814W images that 
Object L could be part of a bar in the lensing galaxy, oriented 
about along the N-S axis. This is consistent with our finding that 
the light profile of the lens galaxy can only be modeled with a 
2-components distribution that includes a significant disk. 

We find no compact and luminous substructure that explains 
the flux ratio anomaly in SDSS J0924-I-0219. However the time 
scale of the anomaly is very large (7 years at least). Therefore 
it seems difficult to explain it only with micro-lensing. A com- 
bination of milli-lensing and micro-lensing is the most likely 
interpretation for the observations. Ob servations of this system 
from UV to Mid-IR and, such as in (lFadelv&Keetonll20Tl1) . 
would probably allow one to conclude about the origin of the 
flux anomaly in SDSS J0924+0219. 
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